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Giotto Ion Mass Spect rometer/High Energy Range Spectrometer (IMS/HERS) 
observations of solar wind ions show charge exchange effects and solar 
wind compositional changes in the coma of comet Halley. As the comet was 
approached, the He to proton density ratio increased until about 1 hour 
before closest approach after which time it decreased. Abrupt increases in 
this ratio were also observed in the beginning and near the end of the so- 
called Mystery Region (8.6 to 5.5x10^ km from the comet along the space- 
craft trajectory) . These abrupt increases in the density ratio were well 
correlated with enhanced fluxes of keV electrons as measured by the Giotto 
plasma electron spectrometer. The general increase and then decrease of 
the He to proton density ratio is quantitatively consistent with a com- 
bination of the addition of protons of cometary origin to the plasma and 
loss of plasma through charge exchange of protons and He^*. In general 
agreement with the solar wind proton and He^ + observations, solar wind 
oxygen and carbon ions were observed to charge exchange from higher to 
lower charge states with decreasing distance to the comet. The more abrupt 
increases in the He 2+ to proton and the He^ + to 0^ + density ratios in the 
mystery region require a change in the solar wind ion composition in this 
region while the correlation with energetic electrons indicates processes 
associated with the comet. 


INTRODUCTION 

Soon after the spacecraft encounters with comet Halley, it was reported that a consid- 
erable amount of solar wind He 2+ was charge exchanged to He + in the comet coma [Balsiger 
etal., 1986; Shelley et al . , 1986]. Specifically, it was reported that inside the mag- 
netic pile-up boundary (1.35 x 10 5 km from the nucleus along the Giotto trajectory) , 
over 30% of the solar wind He^* distribution was charge exchanged to He + . This result 



is illustrated in Figure 1, which shows the time profile of the ratio of M/Q=4 to He 2 * 

+ M/Q=4 densities. This time profile is similar to the one reported by Shelley et al . 
[1986, Figure 3] except that here, the contribution of cometary H 2 + to the M/Q-2 mass 
peak has been removed (see Fuselier et al . [1988] ) . Because of the relatively high con- 
centrations, the M/Q-4 mass peak is thought to contain mainly He* charge exchanged from 
solar wind He 2 * near the comet (within 10^ km from the nucleus) , Figure 1 is expected to 
reasonably represent the fraction of charge exchanged solar wind He 2 * as a function of 
distance from the comet along the Giotto spacecraft trajectory. This figure clearly 
shows that the charge exchanged fraction increased rather abruptly from <10% to over 
40% after the magnetic pile-up boundary (MPB) . (Model results will be discussed in a 
later section . ) 

Solar wind ion charge exchange with cometary neutrals in the comet coma was predicted 
prior to the spacecraft encounters with comet Halley [e.g. , Ip and Ax ford, 1982] . How- 
ever, charge exchange effects were expected to be important at distances less than 
about 10 4 km from the nucleus along the comet-sun line for a Halley-type comet [Ip and 
Hsieh, 1982] . The observations of significant charge exchange of solar wind He 2 * at the 
magnetic pile-up boundary (Figure 1) correspond to a distance that is approximately an 
order of magnitude further from the comet nucleus than predicted. 

Thus far, direct observations of charge exchange of solar wind ions have been limited 
to charge exchange of He^ to He in the comet Halley coma and possibly considerably 
tailward of the nucleus [Milhalov et al , , 1987]. Although the charge exchange cross 
section for He 2 * in H 2 O is not known, it is probably on the order of 3x10“^ cm 2 for 
typical He energies near the magnetic pile-up boundary [Fite et al. , 1962; Koopman, 
1968]. Cross sections for other solar wind minor ions (e.g. O 7 * , O 6 * , etc . ) in H 2 O are 
also poorly known; however, they may be as much as an order of magnitude larger than that 
for He 2 * and nearly equal to the cross section for H* in H 2 0 [ Koopman , 1968; Spjeldvik, 
1979] . Since such a large fraction of He 2 * is charge exchanged in the vicinity of the 
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magnetic pile-up boundary {Figure 1) , it follows that an even larger fraction of other 
solar wind ions should be charge exchanged in that region of the comet Halley coma . The 
larger cross sections of these other solar wind ions should also result in observable 
charge exchange effects further from the comet. 

In this paper, we present observations of solar wind electrons and ions in the coma of 
comet Halley, These observations have potentially important consequences for charge 
exchange near the magnetic pile-up boundary and other locations throughout the comet 
Halley coma . 

Ion observations in this paper were from the Giotto Ion Mass Spectrometer/High Energy 
Range Spectrometer (IMS/HERS) [Balsiger et al . , 1987], This instrument cycled through 
four modes during the Halley encounter. The two modes of interest here measured the 
three dimensional velocity distributions of protons and solar wind ions with mass per 
charge (M/Q)=2 to 4 amu/e. Each individual mode required four seconds (one spacecraft 
spin) to complete an energy and angle scan and the modes were repeated sequentially ev- 
ery 16 seconds. The f ield-of-view of this instrument extended from 15 to 75° relative to 
the spacecraft spin axis (which was also approximately the spacecraft velocity vector 
relative to the comet) and ion energies from 10 eV/e to ^4 keV/e were measured. 

2 i 

For protons and He , the moment calculations have been revised from previous work 
to include an estimate of the percentage of ions outside of the HERS f ield-of-view . 

In this procedure, the distribution was assumed to be gyrotropic in the plasma rest 
frame and the measured phase space density at each velocity and pitch angle was weighted 
by the inverse fraction of phase space sampled at that pitch angle. This procedure 
gave reliable estimates of the density as long as the bulk of the distribution was in 
the f ield-of-view . For protons and He^ + , the bulk of the distribution was in the HERS 
f ield-of-view from upstream of the shock (-10 6 km) to well inside the magnetic pile-up 
boundary (~10 4 km) [Goldstein et al . , 1991], 
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Electron observations in this paper were from the Electron Electrostatic Analyzer 
(RPA1-EESA) [R&me et al . , 1986] . This instrument measured a full 3-dimensional elec- 
tron distribution from 10 eV to 30 keV every 2 seconds. Here, 4 second averages of these 
data were used. 


OBSERVATIONS 

Figure 2 shows the M/Q=2 to proton density ratio at 64 s resolution from 1830 UT to 
near closest approach to the comet. The distance to the comet along the spacecraft 
trajectory is shown above the density ratio profile and some of the major regions and 
boundaries are shown in the figure. The dashed line in this figure will be discussed in 
detail in the model section. 

Upstream from the shock (1830-1922 UT) , the M/Q-2 mass peak is dominated by solar wind 
He and the M/Q-2 to proton density ratio shows rapid fluctuations ranging from 1 . 8 to 
3.8%. These fluctuations are due in part to M/Q-2 counting statistics and the fact that 
the M/Q-2 ions were systematically measured 4 s before the protons. The magnitude of 
the fluctuations gives some indication of the uncertainty in the measured density ratio 
in this region. Downstream from the shock (1922 - 2400 UT) , the average density ratio 
increases from about 2.5% to ~4% before decreasing near closest approach . More abrupt , 
longer lasting deviations are superposed on this trend with the two most notable devia- 
tions occurring at the beginning and near the end of the mystery region. After 2310 UT, 
both the M/Q-2 distribution and the proton distribution are contaminated by a substan- 
tial fraction of cometary ions. The M/Q=2 mass peak contains a significant fraction of 
cometary H 2 + [Fuselier et al., 1988] and the proton peak contains a significant frac- 
tion of cometary pick up protons. While it would be extremely difficult to separate the 
cometary and solar wind protons, the slight difference in the mass per charge of solar 
wind He^ and cometary H 2 has been used to separate these two species [Fuselier et al. , 
1988] . Using the density estimates for solar wind He^ + inside 2310 UT from Fuselier et 
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al . / 5 min averages of the He z to proton density ratio are shown by the filled circles 
connected by the solid line in Figure 2. As can be seen, the He 2+ to proton density ratio 
actually reaches a peak, at about 2310 UT and decreases significantly thereafter. 

In addition to protons and He 2 *, the Giotto IMS/HERS sensor measured solar wind ions 
in the mass per charge range from 2 to 4 . Figure 3 shows two mass spectra from M/Q-2 to 4 
in the solar wind (upper panel) and near the comet (lower panel) . Long averaging times 
were required to obtain good counting statistics for solar wind ions other than He 2 * . 
Plotted is the countrate multiplied by (M/Q) 4 normalized to the peak countrate as a 
function of IMS/HERS mass channel . Error bars are one sigma based on counting statis- 
tics only. The Y-axis is roughly proportional to flux relative to the M/Q-2 ion flux 
(primarily He 2 * in the solar wind and a mixture of He 2 * and H 2 * near the comet) . In the 
solar wind spectrum, the O 7 *, C 5 *, and Ne 8+ fluxes are about 0 . 51 of the solar wind He 2 * 
flux while the O 6 * and C 4 * fluxes are about 1% of He 2 * . In addition, there is a rela- 
tively small amount of M/Q— 4 and slightly larger M/Q ions (probably of solar wind origin 
and consisting mainly of Si 7 * with contributions from multiply charge iron ions) . 

The spectrum in the lower panel of Figure 3 near the comet shows considerable change 
from the solar wind spectrum. Solar wind O 7 *, O 6 *, and Ne 8 * fluxes are all reduced rela- 
tive to the M/Q-2 flux. The most dramatic change occurs in the M/Q-4 flux. It is now well 
above background at a few percent of the M/Q-2 flux. 

Besides M/Q-2, the relatively good counting statistics for O 6 *, the second most abun- 
dant solar wind minor ion, and M/Q-4 ions allow us to construct density profiles for 
these ion species with the time resolution needed to distinguish some features in the 
cometary coma. Figure 4 shows density profiles for these two ion species from the up- 
stream solar wind to near closest approach. A variety of averaging intervals were used 
based on counting statistics of the individual ion species. For example, one hour av- 
erages were used in the solar wind and after 2200 UT and half hour averages were used be- 
tween the shock and the end of the mystery region. Because of the 4 keV/e upper energy 
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per charge cutoff of the IMS/HERS detector, a fraction of the higher mass per charge 
solar wind ion distributions were outside the energy per charge range of the instru- 
ment . This energy per charge cutoff had the greatest effect on the M/Q=4 densities in 
the solar wind, where the solar wind velocity was relatively high. Because mass load- 
ing caused the solar wind velocity to decrease from about 400 km/s far from the comet 
to near zero near the comet, the upper energy per charge cutoff of the HERS detector had 
less effect on M/Q=4 densities closer to the comet. It is estimated that 47% of the dis- 
tribution was above 4 keV/e from 1600-1700 UT, 20% from 1700-1800 UT and about 9% from 
1800-1900 UT . After 1930 UT, the solar wind velocity had decreased enough so that a neg- 
ligible amount of M/Q=4 solar wind ion distribution was above the 4 keV/e cutoff. For 
other solar wind ions, such as C 4+ and 0 5+ , the upper energy per charge cutoff had much 
less effect on the density. It is estimated that less than 5% of the solar wind C 4+ and 
0 5+ distributions were above the 4 keV/e cutoff from 1600-1700 UT . This percentage de- 
creased as the comet was approached. 

Despite the long averaging times, some general trends can be seen in Figure 4 . The 
density profile is similar to the proton and M/Q=2 density profiles [see Goldstein et 
al • / 1987; 1991] . in particular, it shows a factor of two increase in density between 
2000 and 2200 UT associated with the mystery region. In contrast, The M/Q=4 density is 
at detection threshold (~10“^ cm"^) in the solar wind and increases sharply after 2200 
UT [see also, Shelley et al . , 1986]. 

Although the densities of other solar wind ions were extremely low, a general trend 
towards lower charge states was seen in the comet coma. Table 1 lists density ratios for 
Oxygen and Carbon charge states for the 2 hours prior to the crossing of the cometary bow 
shock (^10 6 from the comet) and for the 2 hours before closest approach (~~few x 10^ km 
from the comet) . Whereas the to density ratio decreased with decreasing dis- 
tance to the comet, the to density ratio increased with decreasing distance. 
Also, the C 4+ to C 5+ ratio increased with decreasing distance to the comet. Although 
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the uncertainties are quite large, the general trend in Table 1 was toward lower charge 
states as the comet was approached. 


CHARGE EXCHANGE 

The observations in Figures 1 through 4 and Table 1 are qualitatively consistent with 
charge exchange of solar wind ions in the comet Halley coma. Figure 1 shows that He 2 ' 1 " 
was charge exchanged to He* as the comet was approached. Figure 2 shows a general in- 
crease m the He^ to proton density ratio as the comet was approached. This density ra- 
tio was affected by solar wind composition changes, addition of cometary protons, and 
charge exchange losses. The addition of comeatry protons would cause the He 2 * to pro- 
ton density ratio to decrease with decreasing distance to the comet. It is clear from 
Figure 2 that the presence of cometary protons does not have an effect before 2310 UT, 
when the density ratio began to decrease. The increase in the He 2 * to proton density ra- 
tio from 1830 UT to 2310 UT could be due to a slow change in the solar wind composition 
(i .e. , unrelated to the presence of the comet) . However, the slow change in composition 
would not account for the increase in the M/Q - 4 ion density (interpreted as He*) as the 
comet was approached (Figure 4) . Charge exchange of He 2 * and protons can account for 
the increase m the He^ to proton density ratio in Figure 2 since, as discussed in the 
introduction, the charge exchange cross section for protons is believed to be ten times 
higher than that for He 2 * . Charge exchange of He 2 "*" to He* also accounts for the increase 
in the M/Q=4 ion density as the comet was approached. Although there are several other 
possible sources for the M/Q-4 mass peak (see Figure 3) , the total density of the solar 
wind ions between M/Q-2 to 4 is only a few percent of the He 2 * density. Therefore, these 
ions cannot contribute significantly to the M/Q-4 ion density near the comet and the ma- 
jor contribution to the M/Q-4 mass peak appears to be He* charge exchanged from solar 
wind He 2 * [see also, Shelley et al . , 1986). Superposed on the charge exchange loss of 
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He and protons are possible solar wind composition changes (for example from 2120 to 
2150 UT in Figure 1 ) which will be discussed in more detail in the next two sections . 

Figure 3 and Table 1 are also consistent with charge exchange as the comet was ap- 
proached. In Figure 3, there is evidence for a loss of solar wind minor ions such as 0 7 + 
and Ne° close to the comet. In Table 1, there is a general trend to lower oxygen and 
carbon charge states as the comet was approached. Since all charge states of solar wind 
oxygen have approximately the same charge exchange cross section [Spjeldvi k, 1979 ] , 
the general trend to lower charge states indicates that multiply charged ions undergo 
charge exchange as the solar wind approaches the comet. 

MODEL FOR PROTON AND He 2+ CHARGE EXCHANGE 

Unfortunately, the very low densities of the solar wind minor ions between M/Q=*2 and 
4 allow only a qualitative comparison with predictions from solar wind ion charge ex- 
change. However, the proton, M/Q-2, and M/Q=4 densities are large enough to allow some 
quantitative comparison with predictions from solar wind ion charge exchange. Compar- 
ison between the observations and a simple model for the solar wind interaction have al- 
ready been made [Shelley et al . , 1986]. Shelley etal. concluded that the helium obser- 
vations in the comet coma were inconsistent with present models of the solar wind inter- 
action with the comet. They suggested that model results and observations would be in 
better agreement if the He^ to He charge exchange cross section in H 2 O was substan- 
tially larger than the ~3xl0~ 16 cm 2 assumed in their model and/or if the flow field in 
the comet coma was substantially different than assumed in their model. 

Given a cometary flow model, the ratio of the He 2 * and proton densities can be deter- 
mined along the Giotto trajectory. This ratio will depend on the ratio of the charge ex- 
change cross sections of He^ and protons and on the production rate of cometary H + from 
H 2 O. Since the proton charge exchange cross section in H 2 O at the energies of interest 
here is reasonably well known and the production of H + can be modeled, the profile of the 
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He 2+ to proton density ratio is a function of the He 2 + charge exchange cross section. In 
this section, we use a cometary flow model to predict the He 2+ to proton density ratio 
and compare this prediction to the observations. Through this comparison, we estimate 
the charge exchange cross section for He 2+ in H 2 0. 

The model charge exchange calculations were made using the procedure outlined in Ip 
[1989] . The MHD flow dynamics of the comet-solar wind interaction were separated from 
the photochemistry and the charge exchange process. In doing this, the continuity equa- 
tion 


1 Z ln * v * A) 


( 1 ) 


for the j th species can be integrated. Here, v, is the flow velocity along the stream 
line, n^ is the number density, q j is the production rate, s, is the loss rate, and A is 
the cross section of the stream tube. The flow field model determines the values of v ; 
for each stream line and the cross section of the stream line. Here, we use the Fedder et 
al. 2-D flow model [Fedder et al . , 1986]. 

To compute the net production rate for cometary H+ (q^-s# in Equation 1) , we assume a 
coma model of water vapor and its photodissociation fragments (OH, 0, and H) . The pro- 
duction rate for hydrogen ions is then given by, 


- n{H) {l/ti + y^rijCTj [H) < Vj >) (2) 

i 

Here the ionization time (t*) , which includes photoionization and charge exchange with 
solar wind protons, is assumed to be 10^ seconds , The cross section <Tj (H) in the second 
term refers to the interaction between H + and the neutral gas and <v ; > is the root-mean- 
square of the plasma flow speed v^ and the thermal speed (tfvj) of the j** 1 species . 
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The loss term in (1) is dominated by the charge exchange process. Although charge ex- 
change cross sections for different reactions should vary with the flow speed, in the 
present approximation we assume the the charge exchange rate for protons is constant 
with energy . 

The dashed line in Figure 2 shows the results from the model . In the model, the ini- 
tial density ratio was 2.5%, the He 2 + charge exchange cross section was a = 3x10“^ cm 2 , 
and the H + charge exchange cross section was 5xl0" 15 cm 2 . The general trend of the ob- 
served density ratio is reproduced, indicating that the charge exchange cross section 
for He 2 + is about a factor of 10 lower than that for H + . Increasing the He 2 + charge ex- 
change cross section by a factor of 10 would result in a He 2 + to proton density ratio 
that neither has an initial increase with decreasing distance nor has a decrease in the 
last hour before closest approach. 

Since the densities are computed individually in the model, the He + to total Helium 
ion density ratio can be directly compared with the observations in Figure 1. The per- 
cent charge exchanged He 2 + predicted from the model (solid line connected by x' s in 
Figure 1) is clearly in good agreement with the observations up to the magnetic pile- 
up boundary. This is again consistent with a He^ + charge exchange cross section of 
~3xl0 cm for this region. After the magnetic pile-up boundary , the predicted per- 
centage is about a factor of 4 below the observed percentage. 

Thus, while the general trend in Figures 1 and 2 are reproduced with the expected cross 
sections, the model underestimates the density ratio in the mystery region and near the 
magnetic pile-up boundary, and it also underestimates the percent charge exchanged 
He near the magnetic pile-up boundary. In the next section, we address the discrep- 
ancy between observations and predictions in the mystery region by considering why the 
density ratio changes rather abruptly. The deviation between observations and predic- 
tions is left for the final discussion secion. 
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MYSTERY REGION 


The mystery region is a region of increased density, increased flow speed, decreased 
ion temperature, and increased electron temperature relative to the rest of the re- 
gion downstream from the Giotto bowshock crossing [ Goldstein et al . , 1986; R&me et al . , 
1986; 1990] . Figure 2 also shows that the M/Q—2 to proton density ratio deviates 

from its relatively smooth, increasing trend in this region. In this section, we will 
discuss the possible origins of these deviations and other properties of the mystery 
region. 

Figure 5 shows the density profiles for electrons from 10 eV to 3 . 5 keV, protons, M/Q=2 
ions (primarily solar wind He 2 *) , solar wind O 6 *, and electrons from 0.8 to 3.5 keV. The 
proton and electron density profiles in the top two panels show both qualitative and 
quantitative agreement throughout this time interval. The mystery region extends from 
the sharp density increase at 2022 UT to the sharp decrease in the M/Q-2 ion density and 
the 0.8 to 3.5 keV electron density at 2152 UT. 

The horizontal bars in Figure 5 show the time intervals when the M/Q~2 to proton den- 
sity ratio in Figure 2 deviated from its smooth trend. Of particular interest is the in- 
terval from 2115 to 2150 UT. As can be seen in Figure 5 at 2120 UT, the increase in the 
M/Q-2 to proton density ratio in Figure 2 is associated with a drop in both the elec- 
tron and proton number densities. The M/Q=2 number density did not decrease from 2115 
to 2125 UT, in fact, it increased somewhat. The net result was an abrupt increase in the 
M/Q=2 to proton density ratio. It is also interesting to note that the changes in the 
M/Q=2 to proton density ratio in the period from 2115 to 2150 UT are associated with a 
decrease in the solar wind O 6 * density and an increase in the electron density above 900 
eV (see the bottom two panels in Figure 5) . The half hour averages of the O 6 * density do 
not allow us to determine if a similar decrease in that density is associated with the 
first interval of increased M/Q~2 to proton density ratio from 2013 to 2021 UT . However, 
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it is clear that the increase in the M/Q=2 to proton density ratio is associated with in 
creased electron density above 800 eV from 2037 to 2105 UT. The energetic electron den- 
sity increases in the mystery region are unique . From the bow shock to the mystery re- 
gion and after the mystery region to closest approach, the energetic electron density 
was ^3-4xl0“ 4 cm” 3 (see Figure 5 from 2000 to 2040 UT and R&me et al. [1986] ) . 

DISCUSSION 

In most regions of the comet Halley coma, we have attempted to interpret the IMS/HERS 
observations in terms of charge exchange with cometary neutrals. The increase in the 
He* density with decreasing distance to the comet in Figure 1, the general increasing 
M/Q=2 to H* density ratio in Figure 2, and the trend to lower charge states for multiply 
charged Oxygen and Carbon in Table 1 are all considered as evidence for this charge ex- 
change process . 

While we can find substantial qualitative evidence for charge exchange in the comet 
Halley coma, quantitative comparisons oetween the observations and predictions from 
cometary models is clearly lacking in the mystery region and the region near and inside 
the magnetic pile-up boundary. 

To the properties of the mystery region that have been known previously [e.g. , Gold- 
stein et al. , 1986; R&me et al . , 1990] , we add that there are solar wind composition 
changes in parts of this region. The He 2 * to proton density ratio shows two deviations 
in its increasing trend in this region (see Figure 2) . The second of these deviations 
occurs within the last 30 minutes of the mystery region and is associated with a de- 
crease in the proton density and a possible slight increase in the M/Q*2 density (see 
Figure 5) . Also, it is associated with a decrease in the O 6 * density and the largest 
fluxes of energetic electrons observed during the encounter. 

The general increase in all solar wind ion densities and the change in the He 2f to pro 
ton density ratio are strong indicators that the mystery region is dominated by solar 
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wind plasma with different characteristics from the plasma in other parts of the coma. 
While it is tempting to interpret this region as simply a separate solar wind plasma 
that has convected into the Halley coma, the observation of keV electrons in this re- 
gion indicates that there is also an additional interaction related to the presence of 
the comet. It is also interesting to note that the second increase in the He 2 * to proton 
density ratio (from 2115 to 2150 UT in Figure 5) results primarily from the decrease in 
the proton density and there is a simultaneous decrease in the O 6 * density. As pointed 
out in the introduction, H* and O^* both have charge exchange cross sections that are an 
order of magnitude larger than that for He 2 * . This observation suggests that enhanced 
charge exchange may be occurring in the parts of the mystery region where the He 2 * to 
H + density ratio is elevated and keV electrons are present. The possible causes of this 
enhanced charge exchange and the relation, if any, with the keV electrons are not known 
at this time. Also, it has been noted that the transition from the mystery region back 
into lower density plasma (at 2150 UT in Figure 5) may have been seen by other spacecraft 
that encountered comet Halley and by the Iternational Cometary Explorer (ICE) space- 
craft at comet Giacobini-Zinner [R4>me et al . , 1986; R&me, 1990) . 

Thus, while the observations in Figures 2 and 5 indicate that the mystery region is 
likely a separate solar wind plasma that has convected into comet coma, some of the fea- 
tures in this region and possibly its boundaries indicate that additional and different 
interactions between this plasma and the comet coma are taking place. 

Another important region where the cometary models fail to predict the amount of 
charge exchange is the region from 2310 UT to closest approach, or the region near and 
inside the magnetic pile-up boundary. Figure 1 shows that the amount of charge ex- 
changed He 2 * to He* increases dramatically across the magnetic pile-up boundary. This 

amount is much larger than that predicted by previous models [ Shelley et al . , 1986] and 

by the model used in this paper. In addition, the predicted He 2 * to H* density ratio is 
clearly too low (Figure 2) . 
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The difference between the predicted and observed density ratios indicates that ei- 
ther there is too little He 2 * or too much H* present in the vicinity of the magnetic 
pile-up boundary. As pointed out by Shelley et al. [1986] the sharp increase in the per- 
cent charge exchanged He 2 * at the Magnetic Pile-up Boundary (Figure 1) is not the result 
of a rapid increase in the He* density but a rapid decrease in the He 2 * density. These 
results taken together suggest that the charge exchange of He 2 * relative to H* is under- 
estimated in the cometary models. An overestimate of the the production of cometary H* 
could also help explain the differences between predictions and observations in Figure 
2 in the vicinity of the Magnetic Pile-up Boundary but obviously do not affect the re- 
sults in Figure 1. 

A possibility already discussed by Shelley et al. [1986] is that the He 2 * cross sec- 
tion could be about a factor of 3 to 10 times larger than the 3x10“^ cm 2 value used in 
the models. This is clearly not the case prior to 2310 UT, where the predicted density 
ratio using this cross section agrees well with the observations. One possibility is 
that the cross section for He 2 * increases with decreasing energy faster than that for 
H* . (The H* cross section has been shown to increase with decreasing energy for ener- 
gies below 100 eV [ Koopman 1968] . ) Armed with such a free parameter, the density ra- 
tio profile in Figure 2 could be reproduced exactly, but the understanding of the den- 
sity ratio decrease after 2310 UT would not be improved. Also, this possibility would 
not explain the rather abrupt increase in the percent charge exchanged He 2 * in Figure 
1 because the He * flow velocity decreases smoothly across the magnetic pile-up bound- 
ary and the thermal speed does not change at all across this boundary [Fuselier et al. , 
1987] . One thing that does change rather abruptly in the vicinity of the magnetic pile- 
up boundary is the flow direction [Fuselier et al . , 1987] . 

It is possible that a combination of rapid slowing and deflection of the solar wind 
plasma incident along the sun-comet line, followed by a re-acceleration along the flanks 
of the comet could explain both the large amount of charge exchange and the relatively 
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smooth velocity profile observed by Giotto [ Shelley et al . , 1986; Wallis , 1990}. How- 
ever, the present computer models cannot be used to predict such a possibility. Thus, 
we must conclude that present computer models allow us to qualitatively understand the 
charge exchange of solar wind ions in the outer reaches of the coma and help to distin- 
guish real charge exchange effects from changes in the solar wind composition (for ex- 
ample in the mystery region) . The quantitative understanding of the observations of 
charge exchange in the inner coma (i.e. in the region just before and inside the mag- 
netic pile-up boundary in Figures 1 and 2) remains poor. 
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Table 1. Solar wind Oxygen and Carbon density ratios 


Time (UT) 


0 7 +/ 0 6 + 


0 5 + / 0 6+ 


C 4 +/ C 5 + 


1705-1913 

2200-2400 


0 . 35 ± 0.14 

0 . 17 ± 0.31 


0 . 24 ± 0 . 20 
0 . 59 ± 0 . 81 


0 . 4 2 ± 0 . 4 4 
0 ♦ 74 ± 0 . 92 


If 



Fig. 1. Observed ( + 's) and modeled (x's) M/Q=4 to He 2+ + M/Q=4 density ratio versus dis 
tance from the comet along the Giotto trajectory. The abrupt increase at the magnetic 
pile-up boundary (MPB) is interpreted as a large increase in the percent He 2 + charge ex 
change to He + . 

Fig. 2. M/Q=2 to H + density ratio versus time or distance from the comet {solid line) . 
The filled circles connected by the solid line are the ratios after the cometary H 2 + 
contribution to the M/Q=2 mass peak is removed. The Dashed line show a fit to the ob- 
served ratio from the charge exchange model. 

Fig. 3. Mass spectra from M/Q=2 to 4 far from the comet (upper panel) and near the comet 
(lower panel) . The vertical axis is proportional to the flux normalized by the He 2 + 
flux. Near the comet, the relative fluxes of all multiply charged solar wind species 
are reduced except for M/Q=4 when compared to the mass spectrum further away. 

Fig. 4. Densities versus time for M/Q~4 ions (dashed line) and solar wind 0^ + (solid 
line) . 

Fig. 5. Electron, proton, He 2 "*", 0^"*", and electron >1 keV densities for a time period 
centered on the mystery region. The bars in the figure identify the intervals when the 
He to H density ratio in Figure 2 are elevated. These intervals are correlated with 
enhanced energetic electron flux. The second interval is also correlated with a de- 
crease in the H + and 0 6+ densities . 
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